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y—i ' Abstract. PHENIX measures leptons at mid and forward rapidities and ex- 

. tracts leptons resulting from semi-leptonic decays of heavy quarks. We present 

' PHENIX results related to heavy quark production, specifically the invariant 

' ' . cross section of the non-photonic single electrons produced at midrapidity in 

I p+p and d+Au collisions at ^/snn = 200 GeV . The measured cross section 

■ for p+p collisions is compared with the NLO prediction, and a possible excess 
I/"") ' over the prediction is noted. The measured cross section for d+Au collisions 

scales with the number of colliding nucleon pairs from the p+p spectra. 
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1. Introduction 

•rH , 

■ The parton model and perturbative Quantum ChromoDynamics ( pQCD ) can be 
applied to the violent strong interaction processes characterized by a large scale [ 

As an example, it is applied to the hard particle production where px serves as 
the large scale. Neutral pion production in p+p collisions at y/s = 200 GeV in the 
range 2 < pr < 13 GeV jc measured by the PHENIX experiment at RHIC can be 
well described by the Next-to- Leading-Order (NLO) pQCD [EJ. 

Heavy quark production belongs to the category of violent strong interaction 
processes due to its large mass. While classification of the charm quark as a heavy 
quark is not yet fully settled the parton model and pQCD are often applied to 
the charm quark production. Within the given limitations, predictions for charm 
hadro-production exist from lower collision energies including the Fermi Lab fixed 
target program 

Further predictions were made for heavy quark production in the ion collisions. 
Systematic studies of charm production in p+p and p+ nucleus collisions have been 
proposed as a sensitive way to measure the parton distribution function in nucleons, 



2 



Y. Kwon et al. 



and nuclear shadowing effects [0]. For heavy ion collisions at RHIC energies, heavy 
quark energy loss [ 5 , modification of charm quark hadronization [ |B] , in-medium 
effects on open- and hidden-charm states [0, possible J/ip suppression and 
charm flow have been proposed as possibilities. 

2. Theoretical background 

Factorization of the hadro-production cross section of heavy flavored hadrons is 
described as follows. The differential cross-section for the production of a heavy 
flavored particle (a charmed particle in the particular case), H, in the collision of 
two hadrons, A and B, is believed to factorize in the following manner 

da[A + B -> H + X] = ^ ft/ A® fj/B®da[ij -> hh + X](g> D h ^ H + remainder (1) 

Here i,j denote partons which are point-like, and fj/ B are parton distri- 
bution functions, and Dh^H is the fragmentation function for the heavy quark h 
hadronizing into H. The parton distribution functions and the fragmentation func- 
tions are process independent. The da[ij — > hh + X] is a perturbatively calculable 
short-distance cross section. The remainder represents corrections to the factor- 
ized form of the cross section that are power-suppressed by AQCD/rrih or possibly 
A-QCd/pt if Pt ^> mil- A large quark mass nih makes the remainder smaller, and 
application of pQCD to the factorized term becomes possible. 

As the collision energy increases, the factorized term increases faster than the 
remainder, and the hard processes (factorized term) gain importance at high colli- 
sion energies. We can test the validity of the factorization at RHIC collision energies 
by measuring heavy flavored hadron production for p+p collisions. 

The factorization scheme has further consequences for the collisions of ions. For 
the interaction between point-like partons, we can approximate the nuclear parton 
distribution as follows. 

fi/Au ~ 79 fi/ p + 118 /,/„ ps 197 /j/jv, 

and 

fi/d ~ fi/p + fi/n ~ 2 fi/N 

where /j/at is the parton distribution function inside the nucleon. 

If the charm producing processes are point-like and there is no modification 
of the initial parton distribution for the d + Au collisions, the charm production 
cross section will scale with 2 x 197, i.e. the number of colliding nucleon pairs. 
Surprisingly, this extension does not work for high pt particles produced in central 
AU+ Au collisions [ ^] , and careful studies of the stated scaling are desirable for 
the various hard processes. 
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3. Measurements 



PHENIX is capable of measuring leptons at mid and forward rapidities. PHENIX 
uses global detectors to characterize the collisions, a pair of central spectrometers 
at mid rapidity to measure electrons, hadrons, and photons, and a pair of forward 
spectrometers to measure-^riuaiis. Each spectrometer has a large geometric ac- 
ceptance of about one fpliJiilipMiiyt'ellent energy ^ijd^g^OQ^p^^ni resolution, and 
particle identification . (see ^ T ' ' l f " 
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Fig. 1. Acceptance of the PHENIX experiment. 



The PHENIX central arms consists of tracking systems for charged particles 
and electromagnetic calorimetry. The calorimeter is the outermost subsystem of 
the central arms and provides measurements of both photons and energetic elec- 
trons. The tracking system uses three sets of Pad Chambers (PC) to provide precise 
three-dimensional space points needed for pattern recognition. The precise projec- 
tive tracking of the Drift Chamber (DC) is the basis of the excellent momentum 
resolution. A Time Expansion Chamber (TEC) in the east arm provides additional 
tracking and particle identification. The Time-of-Flight (ToF) and Ring-Imaging 
CHerenkov (RICH) detectors provide particle identification. The RICH provides 
separation of electrons from the large number of copiously produced pions. Using 
information from the RICH, the TEC, and the electromagnetic calorimeter it is 
possible to reject pion contamination of identified electrons to one part in 10 4 over 
a wide range of momenta. 

Each forward spectrometer is based on a muon tracker inside a radial magnetic 
field followed by a muon identifier, both with full azimuthal acceptance. The muon 
trackers consist of three stations of multiplane drift chambers that provide precision 
tracking. The muon identifier consist of alternating layers of steel absorbers and 
low resolution tracking layers of streamer tubes of the Iarocci type. With this com- 
bination, the pion contamination of identified muons is typically 3 x 10~ 3 , and can 
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be estimated and subtracted statistically. Measurement of prompt muons resulting 
from semi-leptonic decays of heavy flavored hadrons is also possible, and analysis is 
actively persued. 
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Fig. 2. The inclusive single electron spectra and composition of backgrounds, p+p 
collisions at ^snn = 200 GeV. 



PHENIX measures the non-photonic single electron spectra by first measuring 
the inclusive electron spectra and then subtracting the contribution from photonic 
sources. Inclusive electrons contain two components: (1) non-photonic - primarily 
semi-leptonic decays of mesons containing heavy (charm and bottom) quarks, and 
(2) photonic - Dalitz decays of light neutral mesons (tt , rj, 77', p, u>, and tfi) and photon 
conversions in the detector material. The PHENIX experiment uses two approaches, 
converter subtraction and cocktail subtraction, to estimate the photonic component. 

The converter subtraction is a data-driven approach. A photon converter (a 
thin brass tube of 1.7% radiation length thick) is installed in the middle of data- 



Heavy flavor production at ^/snn = 200 GeV 



5 



taking. The photon converter multiplies the photonic contribution to the electron 
yield in a well-defined manner and hence the photonic component can be estimated. 

The cocktail subtraction simulates electron production from the known sources 
and subtract them from the inclusive spectra. PHENIX has measured the pr dis- 
tributions of 7r ± and 7T°. We fit a power law function to the combined data sets 
to determine the input 7r° spectrum for the decay generator. The pt distribu- 
tion of any other hadrons is obtained from the 7r° spectrum by replacing pt with 



+ m\ — m^ . In this parameterization h/ir° ratios approach constants at high 
Pt- The asymptotic ratios used to fix the relative normalizations are determined 
from the world data including the PHENIX results. Photon conversions are evalu- 
ated using a combination of the GEANT simulation and the hadron decay generator. 

The converter subtraction and the cocktail subtraction yield consistent results, 
and Fig. [21 shows the inclusive electron spectra and composition of the estimated 
background. 

4. Result 

The cocktail subtracted invariant electron cross section is shown on Fig. [3] The 
PYTHIA calculation (m c = 1.25 GeV/c 2 , m b = 4.1 GeV/c 2 , K = 3.5, and < 
k t > = 1.5)) is displayed on the same figure. The NLO pQCD prediction of heavy 
flavor-related electrons was made using the HVQLIB package. The heavy quark 
production cross section was calculated using HVQLIB and PYTHIA, and the ratio 
of heavy quark yields {dN q / dpx) hv qli b to (dN q /dpr) pythia was used as a weight 
for the PYTHIA decay electrons. This method enables us to use an exact NLO 
pair production calculation together with fragmentation and decay kinematics from 
PYTHIA. Within the NLO pQCD calculation, charm decays are the major source 
of non-photonic single electrons for most of the observed pt range. The NLO 
prediction for open charm yields harder spectra than the PYTHIA calculation, but 
an excess of data over both predictions is seen beyond intermediate px- An elaborate 
theoretical study has appeared recently and also shows the observed excess [ lllj. 
The uncertainty of theoretical prediction mostly comes from the mass and the scale 
uncertainties, and is relatively large. Reduced experimental errors will contribute 
to a definitive statement. 

Fig. 01 shows both the non-photonic single electron production in the d+Au 
collisions at ^/snn — 200 GeV scaled with the number of colliding nucleon pairs, 
described in section|3 and the empirical parametrization oip+p production. We ob- 
serve that the non-photonic single electron production in the d+Au collisions scales 
with the number of colliding nucleon pairs. Further study was made for four cen- 
trality classes. Non-photonic single electron spectra obtained from each centrality 
class are also consistent with the p+p spectra scaled by the nuclear thickness. The 
observation is consistent with the picture of non-photonic single electrons produced 
by point-like interactions. 
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Fig. 3. The non-photonic single electron spectra from the p+p collisions at 
Y / J/vjv = 200 GeV and comparison with the NLO pQCD prediction. 

5. Conclusion 

We have discussed exciting theoretical and experimental aspects of heavy quark 
production with emphasis on quantities measurable by PHENIX. Selected results 
from PHENIX were presented and we note key experimental observations. 

a) The non-photonic single electron spectra from p+p collisions at ^snn = 
200 GeV indicate a possible excess over the NLO pQCD heavy quark prediction. 

b) Within errors, the non-photonic single electron spectra obtained from d+Au 
collisions both in minimum bias or in various centrality classes exhibits scaling with 
the number of colliding nucleon pairs from the p+p spectra. Hence, production 
process for non-photonic single electrons is consistent with point-like interactions. 
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Fig. 4. The non-photonic single electron spectra from the d+Au collisions at 
y/sJiN = 200 GeV (scaled with the number of colliding nucleon pairs). The empir- 
ical parametrization of the p+p spectra is also displayed. 
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